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ABSTRACT 

Context. Diffuse intracluster light (ICL) and cD galaxy halos are believed to originate from galaxy evolution and disruption in 
clusters. 

Aims. The processes involved may be constrained by studying the dynamical state of the ICL and the galaxies in the cluster core. 
Here we present a kinematic study of diffuse light in the Hydra I (Abell 1060) cluster core, using planetary nebulas (PNs) as tracers. 
Methods. We used multi-slit imaging spectroscopy with F0RS2 on VLT-UTl to detect 56 PNs associated with diffuse light in the 
central 100 x 100 kpc^ of the Hydra I cluster, at a distance of ^ 50 Mpc. We measured their [OIII] msoo? magnitudes, sky positions, 
and line-of-sight velocity distribution (LOSVD), and compared with the phase-space distribution of nearby galaxies. 
Results. The luminosity function of the detected PNs is consistent with that expected at a distance of ~ 50 Mpc. Their number density 
is ~ 4 times lower for the light seen than expected, and we discuss ram pressure stripping of the PNs by the hot intracluster medium 
as one of the possible explanations. The LOSVD histogram of the PNs is highly non-Gaussian and multipeaked: it is dominated by a 
broad central component with a ~ 500 km s~^ at around the average velocity of the cluster, and shows two additional narrower peaks 
at 1800 km s^^ and 5000 km s^^ . The main component is broadly consistent with the outward continuation of the intracluster halo of 
NGC 331 1, which was earlier shown to have a velocity dispersion of ~470km s^^ at radii of > 50". Galaxies with velocities in this 
range are absent in the central 100 x 100 kpc'^ and may have been disrupted earlier to build this component. The PNs in the second 
peak in the LOSVD at 5000 km s^^ are coincident spatially and in velocities with a group of dwarf galaxies in the MSIS field. They 
may trace the debris from the ongoing tidal disruption of these galaxies. 

Conclusions. Most of the diffuse light in the core of Abell 1060 is still not phase-mixed. The build-up of ICL and the dynamically 
hot cD halo around NGC 3311 are ongoing, through the accretion of material from galaxies falling into the cluster core and tidally 
interacting with its potential well. 

Key words. galaxies:clusters:general - galaxies:clusters:individual (Hydra I) - galaxiesxD - galaxies:individual (NGC 3311) - 
planetary nebulae: general 

1 . Introduction proach to deep photometry for studying the ICL, also enabling 

,. , ^ . r- , , , kinematic measurements for this very low surface brightness 

Intracluster light (ICL) consists or stars that hll up the cluster noDulation 

space among galaxies and that are not physically bound to any a • . . . • w i. . 

^, , , 1 1 1 • 1 An important open question is the relation between the 

galaxy cluster members. For clusters in the nearby universe, the ^/-^t j » j j 1 u i f u • . i . i • 

^ : , , ... ^ , , , ICL and the extended outer halos or brightest cluster galaxies 

morphology and quantitative photometry or the ICL have been /r>/-^/-^ \ u ■ j j . . u »u 

1 • 1 • 1 1 . 1 11 • ^ • 1 (BCGs), whether they are independent components or whether 

studied with deep photometric data or by detection or single stars r • j- i . • r.u i „ n • i c 

. . j,^/; ■' the former is a radial extension or the latter Using a sample or 

in large areas ot sky. , , . . 683 SDSS clusters, Zibetti et al. (2005) found a surface bright- 

Deep large-neld photometry shows that ICL is common in i /4 

S clusters of galaxies and it has morphological structures with "^^^ ^'^'^^^s '"^^P'^'^t t° ^" i""^'" ^ P^'^^l^ "^^'1 1° '^^^'^"be 

■ ■ different angular- scales. The fraction of light in the ICL with *e mean profile of the BCGs, but it is not known yet whether 

respect to the total light in galaxies is between 10% and ^"^^1°P^ ^^P^^ '^^"^^^ P^"^ °f *^ '^1^^^^'" ^ ^^^f" 

30%, depending on th e cluster mas s and evolutionary status ^""'^ component o r whether it is distinc t from the ICL and 

(■Feldmeieret al' OT; 'Adami et al.' OT; 'Mihos et al.' OT, P^'^ *^ ^ost galaxy ^Gonzalez et al. || 2005|) . 

[Zibetti et al. 2005; Krick & Bernstein 2007; Pierini et al. 2008). ^oth the ICL and the halos of BCGs are believed to have 

The detection of individual stars associated with th e ICL, such fo™ed from stars that were tidally dissolved from their former 

as planetary nebulas (PNs) (Arnaboldi et al. 2004; Ague rrietalJ host galaxies or from entirely disrupted galaxies. A number of 

120051: iGerhard et al.ll2007t fC astro-Rodriguez et al. 2009), glob- processes havebeen discussed, startmg withearly work such as 



ular clusters (GCs) (Hilker 2002; Lee et al. 2010), red giants | Richstone| ([i97i;| Hausman & Ostrikeij ([i97i). Contributions to 

stars (Durrell et al. 2002; WiUiamset al. 2007), and supernovae ICL are thought to come from weakly bound stars gener- 

dGal-Yamet al.ll2503t lNeiIlttlLll20d5h is a complementary ap- ^t^d by interactions in galaxy groups, subsequendy released in 

— the cluster's tidal field (Rudick et al. 2006. 2009: Kapfe rer et all 

Send offprint requests to: G. Ventimiglia, e-mail: gven- 2010), interacti ons of galaxies with each other and with the clus- 

timi@mpe.mpg.de ter's tidal field dMoore et alJll996l:rGnedinll2003l:IWiIlman et al.l 
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12004 . and from tidal dissolution of sta rs from massive gal- 
axies prior to mergers with the BCG (iMurante et alJ I2007t 
[Puchwein et al. 2010). Stars in BCG halos may have originated 
in both such major mergers as well as through minor mergers 
with the BCG. Which of these processes are most important is 
still an open issue. 

Kinematic studies of the ICL and the cD halos are instru- 
mental in answering these questions. The kinematics of the ICL 
contains the fossil records of past interactions, due to the long 
dynamical timescale, and thus helps in reconstructing the pro- 
cesses that dominate the evolution of galax ies in clusters a nd the 
formation of the ICL (Rudick et al. 2006; Gerhard et al. 2007t 
IMurante et al..,2007.; . Arnaboldi & Gerhard 2010.). The kinemat- 
ics in the cD halos can be used to separat e cluster from gal - 
axy components, as shown in simulations dPolae et al.ll2010l) : 
so far, however, the observation al results are not un animous: in 
both NGC 6166 in Abell 2199 (iKelson et alJl2 002') as well as 
NGC 331 1 in Abell 1060 (I Ventimiglia et al.ll2oiota) the velocity 
dispersion profile in the outer halo rises to nearl j>^ cluster val- 
ues, whereas in the Fornax cD galaxy NGC 139 9 ('McNeil et alj 
[2010) and in the central Coma BCGs dCocca to et al. 2010) the 
velocity dispersion profiles remain flat, and in M87 in Virgo 
(iDohertv et al.ll2009i) it appears to fall steeply to the outer edge. 
Evidently, more work is needed both to enlarge the sample and 
to link the results to the evolutionary state of the host clusters. 

The aim of this work is to further study the NGC 3311 
halo, how it blends into the ICL, and what is its dynamical 
status. NGC 3311 is the cD galaxy in the core of the Hydra 
I (Abell 1060) cluster Based on X-ray ev idence, the Hydra I 
cluster is the p rototype of a relaxed cl uster (Tamura et al. 2000; 
iFurusho et aTji2001:.Christlein & Zabludoff 2003) . Surface pho- 
tometry is available in the Johnson B, Gunn g and r bands 
(IVasterberg et al.iri99lh . and the velocity dispersion profile has 
been measured out to 100" (jVentimiglia et al. 2010b), show- 
ing a steep rise to ^ 470 km s^^ in the outer halo. Here we use 
the kinematics of PNs from a region of 100 x 100 kpc^ cen- 
tered on NGC 331 1, to extend the kinematic study to larger radii 
and characterize the dynamical state of the outer halo and of the 
cluster core. 

In Section |2] we summarize the properties of the Hydra I 
cluster from X-ray and optical observations. In Section |3] we 
discuss PNs as kinematical and distance probes, and the "Multi- 
Slit Imaging Spectroscopy - MSIS" technique for their detec- 
tion in clusters in the distance range 40 — lOOMpc. We present 
the observations, data reduction, identification, and photometry 
in Sections |4] and |5] In Section |6] we describe the spatial distri- 
bution, line-of-sight (LOS) velocity distribution (LOSVD), and 
magnitude-velocity plane of the PN sample. In Section |7] we 
use the properties of the planetary nebulae luminosity function 
(PNLF) and a kinematic model for the PN population to pre- 
dict its LOSVD in MSIS observations. The simulation allows 
us to interpret the observed LOSVD and also to determine the 
luminosity-specific PN number or a parameter for the halo of 
NGC 3311. In Section [8] we correlate the velocity subcompo- 
nents in the PN LOSVD with kinematic substructures in the 
Hydra I galaxy distribution and discuss implications for galaxy 
evolution and disruption in the cluster core. Finally, Section |9] 
contains a summary and the conclusions of this work. 



2. The Hydra I cluster of galaxies (Abell 1060) 

The Hydra I cluster (Abell 1060) is an X-ray bright, non- 
cooling flow, medium compact cluster in the southern hemi- 



sphere, whose central region is dominated by a pair of non- 
interacting giant efliptical galaxies, NGC 3311 and NGC 3309. 
NGC 3309 is a regular giant elH ptical (E3) and NGC 33 11 is a 
cD galaxy with an extended halo dVasterberg et al.l[l99lh . 

X-ray properties of Hydra I - Except for two peaks associated 
with the bright efliptical galaxies NGC 3311 and NGC 3309, 
the X-ray emission from the hot intracluster medium (ICM) 
in the Hydra I (A 1060) cluster is smooth and lacks promi- 
nent spatial substructures. The center of the nearly circularly 
symmetric emissio n contours roughly coincides with the cen- 
ter of NGC 3 311 (iTamur a et al.l l200flt lYamasaki et al.1 llool 
iHavakawa et a l. 2004, 2006). A faint extended emission with an- 
gular scale < 1' trailing NGC 331 1 to the northeastern, overlap- 
ping with an Fe excess, could be due to gas stripped from NGC 
3311 if the galaxy moved towards the south-west with veloc- 
ity > 500kms"\ according to ' Havakawa et all ( |2004 l2006l) . 
The total gas mass and iron mass conta ined in this region are 
~ 10^ Mq and 2 x 10^ Mq, respectively (Havakawa et al. '20041 
I2006h . The emission components of NGC 331 1 and NGC 3309 
themselves are small, extending to only ^ 10" ~ 2.5 kpc, sug- 
gesting that both galaxies lost most of their gas in earlier interac- 
tions with the ICM. In both galaxies, the X-ray gas is hotter than 
the equivalent temperature corresponding to the central stellar 
velocity dispe rsions, and in approx imate pressure equilibrium 
with the ICM (lYamasaki et al.ll2002l) . 

On cluster scales the X-ray observations show that the 
hot ICM has a fairly uniform temperature distribution, rang- 
ing from about 3.4 KeV in the center to 2.2 KeV in the outer 
region, and constant metal abundances out to a radius of 230 
kpc. Deviations from uniformity of the hot gas temperature 
and metallicity distribution in Hydra I are in the high metal- 
licity region at ^ 1.5arcmin northeastern of NGC 3311, and 
a region at a slightly higher temperature at 7 arcmin south- 
east of NGC 3 311 (Tamura etal. 2000; Furusho et al. 200ll 
l Yamasaki et"al] 12002; .Havakawa et al.. .2004 ,2006: Sato et al] 
2007|). Based on the overall regular X-ray emission and tem- 
perature profile, the Hydra I cluster is considered as the pro- 
totype of an evolved and dynamically relaxed cluster, with the 
time elapsed since the last major subcluster merger being at 
least several Gyr. From the X-ray data the central distribution 
of dark matter in the cluster has been estimated, giving a cen- 
tral density slope of_^_ — 1.5 and a mass within 100 kpc of 
~ 10" Mq (Tam ura et a U lToOO; Hayaka\va et al. 2004). Given 
these properties, the Hydra I cluster is an interesting target for 
studying the connection between the ICL and the extended halo 
of NGC 3311. 

The cluster average velocity and velocity dispersion - From 
a deep spectroscopic sample of cluster galaxies extending to 
Mr < -14, Christlein & Zabludoff (2003) derive the average 
cluster redshift (mean velocity) and velocity dispersion of Hydra 
I. We adopt their values here: vny = 3683 ± 46kms~^, and 
(Tjjy — 724 ± 31kms^^. The sample of measured galaxy spec- 
tra in Hydra I is extended to fainter magnitudes My > —17 
thro ugh the c a talog of early-type dwarf galaxies published by 
Mis geld et al.l (l2008l) : their values for the average cluster ve- 
locity and velocity dispersion are WHy = 3982 ± 148kms~^ 
and (THy = 784km s~^, with the average cluster velocity at 
somewhat higher valu e with respect to the measurement by 
IChristlein & Zabludof j d2003l) . Both catalogs cover a radial 
range of ~ 300 kpc around NGC 3311. Close to NGC 3311, 
a predominance of velocities redshifted with respect to vny is 
seen, but in the radial range ^ 50 — 300 kpc, the velocity distri- 
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bution appears well-mixed with about constant velocity disper- 
sion. 

Distance estimates - The distance to the Hydra I cluster 
is not well constrained yet, as different techniques provide 
rather different estimates. The cosmological distance to Abell 
1060 based on the cluster reds hift is 51.2 zfc 5.7 Mpc assum - 
ing iJo = 72 ± 8km- 1 Mpc^i dChristlein & Z abludoff 2003), 
while direct measurements using the surface brightness fluctu- 
ation (SBF) method for 16 gala xies give a distance of 41 Mpc 
dMieske. Hilker. & Infantgi2005l) . 

The relative distance of NGC 3311 and NGC 3309 along the 
line of sight is also controversial. Distance measurements based 
on the globular cluster luminosity function locate NGC 3311 
about 10 M pc in front of NGC 3309, which puts NGC 3309 
at 61 Mpc (Hilker"2003'), while SBF measurements suggest the 
opposite, with NGC 3311 now at shorter distance of about 
41 Mpc and N GC 3309 even clo ser at 36 Mpc, 5 Mpc in front 
of NGC 3311 (lMieskeetal.ll2005h . 

In this work we assume a distance for NGC 3311 and the 
Hydra I cluster of 51 Mpc, corresponding to a distance modu- 
lus of 33.54. Then 1" corresponds to 0.247 kpc. The systemic 
velocity for NGC 3311 and its central velocity dispersion are 
VN3311 = 38 2 5 (38 00) ±8kms-^ (heliocentric; without and in 
brackets with relativistic correction), and cto = 154± 16kms-^ 
(fV^timiglia et al.ll2010bh . The systemic v elocity of NGC 3309 
is VN3309 = 40 99 km s"^ (iMisgeldet al.l[2008i) . The velocities 
of the other Hydra I gala xies are extracted from t he cat alogs of 
iMisgeld et al] (l2008h and lChristlein & Zabludofj (l2003h . 



3. Probing the ICL kinematics using pianetary 
nebuias 

3. 1. Planetary nebulas as kinematical probes and 
distance Indicators 

PNs occur as a brief phase during the late evolution of solar-type 
stars. In stellar populations older than 2 Gyrs, about one star 
every few million is expe cted to be in the PN phase at any one 
time dBuzzoni et a l. ■2006^. Stars in the PN phase can be detected 
via their bright emission in the optical [OIII]A 5007 A emission 
line, because the nebular shell re-emits ~ 10% of the UV pho- 
tons e mitted by the stellar core in this single line dCiarduUo et alj 
l2005h . When the [OIII] emission line is detected, the line-of- 
sight velocity of the PN can be easily measured. 

The number density of PNs traces the luminosity density of 
the parent stellar population. According to single stellar popula- 
tion theory, the luminosity-specific stellar death rate is indepen- 
dent of the precise star formation history of the associated stellar 
population dRenzini & B uzzoni 1986; Buz zoni et a lj2006h . This 
property is captured in a simple relation such that 



The planetary nebula luminosity function (PNLF) technique 
is one of the simplest methods for determining extragalactic dis- 
tances. This is based on the observed shape of the PNLF. At faint 
magnitudes, the PNLF has the power-law form predicted from 



NpN = OiLgal 



(1) 



where NpN is the number of all PNs in a stellar population, Lgai 
is the bolometric luminosity of that parent stellar population and 
a is the luminosity-specific PN number. The predictions from 
stellar evolution theory are further supported by empirical evi- 
dence that the PN numbe r density profiles follow light in late- 
and e arly-type galaxies dHerrmann et alj 120081 ; ICoccato et alj 
1200 9^. and that the luminosity-specific PN number a stays more 
or less constant with (B-V) color. The empirical result that the 
rms scatter of a for a gi ven color is about a f actor 2-3 remains to 
be explained, however dBuzzoni et alJl2066l) . 



models of uniformly e xpanding shells sur r oundi n g slowly evolv - 
ing central stars CHe nize & Westerlundl 119631: Ijacobvl [T980h . 



However, observations and simulations have demonstrated that 
the bright end of the PNLF dramatically breaks fr om this relation 
and falls to zero very quickly, w ithin ~ 0.7 mag dCiardullo et al.l 
ll998l;lMendez & Soffneilll997h . It is the constancy of the cutoff 
magnitude, M* = —4.51, and the high monochromatic luminos- 
ity of PNs, that makes the PNLF such a useful standard candle. 

3.2. The Multl-Sllt Imaging Spectroscopy technique 

At the distance of the Hydra I cluster, the brightest PNs at 
the PNLF cutoff have an apparent msoo? magnitude equal to 
29.0, corresponding to a flux in the [OIII]A5007A line of 
8 X 10~^^ erg s~^cm~^ according to the definition of 7715007 by 
Jacobv (1989). To detect these faint emissions we need a tech- 
nique that substantially reduces the noise from the night sky. 
This is possible by using a dedicated spectrosc opic technique 
named "Multi-Slit Imaging Spectroscopy" (MSIS. IOerhard et aP 
2005; Arnaboldi et al. 2007). 

MSIS is a blind search technique that combines the use of a 
mask of parallel slits, a dispersing element, and a narrow band 
filter centered at the redshifted [OIII]A5007A emission line. 
With MSIS exposures, PNs and other emission objects in the 
filter's wavelength range which happen to lie behind the slits are 
detected, and their velocities, positions, and magnitudes can be 
measured at the sam e time. The [ OIII] emission line from a PN 
is ^30kms-^ wide dArnaboldi e t al. 2008), so if dispersed with 
a spectral resolution R ^ 6000, it falls on a small number of 
pixels, depending on the slit width and seeing. 

In this work we use MSIS to locate a sample of PNs in the 
core of the Hydra I cluster and measure their velocities and mag- 
nitudes. Our aim is to infer the dynamical state of the diffuse 
light in the cluster core, as described below in Sections|2]and[8] 



4. Observations 

MSIS data for Hydra I were acquired during the nights of March 
26-28, 2006, with FORS2 on UTl, in visitor mode. The FORS2 
field-of-view (FoV) is ^ 6.8 x 6.8arcmin^, corresponding to 
~ 100 X 100 kpc^ at the distance of the Hydra I cluster The ef- 
fective field area in which it was possible to position slits with 
the Grism used here is 44.6 arcmin'^. The FoV was centered on 
NGC 3311 at a = 10h36m42.8s, S = -27d31m42s (J2000) in 
the core of the cluster. The FoV is imaged onto two 2x2 re- 
binned CCDs, with spatial resolution 0".252 per rebinned-pixel. 
The mask used has 24 x 21 slits, each 0".8 wide and 17". 5 long. 
The area covered with the mask is about 7056arcsec^, corre- 
sponding to about 4.4 % of the effective FoV. To cover as much 
of the field as possible, the mask was stepped 15 times so as to 
fill the distance between two adjacent slits in the mask. The total 
surveyed area is therefore 29.4arcmin^, i.e., 66 % of the effec- 
tive FoV. Three exposures of 800 sec were taken at each mask 
position to facilitate the removal of cosmic rays during the data 
reduction process. 

The dispersing element was GRJSM-600B with a spec- 
tral resolution of 0.75 ApixeP^ (or 1.5 Arebinned-pixeP^) at 
5075 A. With the adopted slit width, the measured spectral res- 
olution is 4.5A or 270 kms^^. Two narrow band filters were 
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used, centered at 5045 A and 5095 A, respectively, both with 
60 A FWHM. This ensures the full coverage of the Hydra I clus- 
ter LOS velocity range. Each illuminated slit in the mask pro- 
duces a two-dimensional spectrum of 40 rebinned pixels in the 
spectral direction and 70 rebinned pixels in the spatial direction. 

The seeing during the observing nights was in the range from 
0".6 to 1".5. For the average seeing (0".9) and with the spectral 
resolution of the set-up, monochromatic point-like sources ap- 
pear in the final spectra as sources with a total width of ~ 5 
pixels in both the spatial and wavelength directions. 

Biases and through-mask flat field images were also taken. 
Arc-lamp calibration frames with mask, Grism and narrow band 
filter were acquired for the extraction of the 2D spectra, their 
wavelength calibration and distortion correction. Long slit data 
for the standard star LTT 7379 with narrow band filter and Grism 
were acquired for flux calibration. 

5. Data reduction and analysis 

The data reduction is ca rried out in IRAF as described in 
lArnaboldi etaTI (l2007h and lVentimi^fia et al.l(l2008l) . The frames 
are registered and co-added after bias subtraction. The con- 
tinuum light from the bright galaxy halos is subtracted us- 
ing a median filtering technique implemented in the IRAF 
task .images. im filter. median, with a rectangular window of 
19 X 35 pixels. Then emission line objects are identified, and 
2D-spectra around the emission line positions are extracted, 
rectified, wavelength and flux calibrated, and background sub- 
tracted. Finally the wavelength of the redshifted [OIII]A 5007 A 
emission line for all the identified sources is measured via a 
Gaussian fit. The heliocentric correction for the PN velocities 
is -5.44 kms"^ 

5. 1. Identification of Emission-Line Objects 

All emission line objects found are classified according to the 
following criteria as 

- PN candidates: unresolved emission line objects, both in 
wavelength and spatial direction, with no continuum; 

- background galaxy candidates: unresolved emission line ob- 
jects with continuum or resolved emission line objects both 
with and without continuum. 

The total number of detected emission line sources in our data 
set is 82, of which 56 are classified as PN candidates and 26 as 
background galaxy candidates, of which 6 are classified as [Oil] 
emitters and the remaining 20 as candidate Lya galaxiefl 

For details on the background galaxy candidates see 
IVent imiglia e t alJ(l2010ah . Note that the background galaxy clas- 
sification is independent of luminosity and that these objects 
have a broad equivalent width distribution. Therefore, the fact 
that the PN candidates (unresolved emission line objects without 
detectable continuum) have a luminosity function as expected 
for PNs observed with MSIS at a distance of ^ 50 Mpc (see 
Section ITTT i. implies that the large majority o f these PN candi- 
dates m ust indeed be PNs. In addition. Fig. I of lVentimiglia et al] 
(1201 Oah shows that all of the background galaxy candidates 

' Note that the equivalent widths (EWs) of the PN candidates are 
mostly distributed between 30 A < EW < 100 A , similar to the EWs 
of the background galaxy candidates, and cannot therefore be used to 
discriminate between both types of emission sources. This is because 
these distant PNs are faint and the continuum level in the MSIS images 
is given by the Icr limit from the sky noise; see Section [5!2l 



but two fall in the blue filter in the velocity range between 
1000 kms"^ and 280Gkms"\ blue-shifted with respect to the 
Hydra I cluster. Several of the unresolved background galaxies 
in this blue-shifted velocity range have a continuum level just 
above the detectability threshold, suggesting that the PN candi- 
date sample may contain a few background galaxy contaminants 
in this velocity range whose continuum is too faint to detect. 

The two background galaxies seen in the red filter are both 
extended and have medium bright emission fluxes; one has a 
very bright continuum, the other no detectable continuum. From 
this we conclude that the residual contamination of the PN candi- 
date sample at velocities > 3000 kms^^ must be minimal. With 
this in mind, we will in the following simply refer to the PN 
candidates as PNs. 

5.2. Pfiotometry 

Magnitudes of the PN candidates are computed using the msoor 
definition by Jacoby (1989), TO5007 = -2.51ogF5oo7 - 13.74, 
where -F5007 is the integrated flux in the line computed in circular 
apertures of radius 0".65 — 0".85 in the 2D spectra, measured us- 
ing the IRAF task .noao.digiphot.aphot.phot. The Icr limit on 
the continuum in these spectra is 7.2 x 10~^*'ergcm^^s^^A 

5.2.1. Photometric errors and completeness function 

The photometric errors are estimated using simulations on a 
sample of 2D wavelength, flux calibrated and background sub- 
tracted spectra. For each simulation 100 artificial PN sources are 
generated using the IRAF task .noao.artdata.mkobject. The 
adopted PSF is a Gaussian with a dispersion obtained by fit- 
ting a 2D Gaussian to the profile of a detected PN candidate 
with adequate signal-to-noise. The a value is 1.1 pixels, i.e., 
FWHM - 0".7, and the FWHM in wavelength is - 4 A. The 
simulated PN samples have luminosity functions (LFs) given by 
a delta function at one of five different input magnitudes (29.3, 
29.7, 30.1, 30.5 and 30.9 mag). The output magnitudes on the 
2D spectra are measured with .noao.digiphot.aphot.phot us- 
ing circular apertures, in the same way as for real PN candidates. 
In these experiments, no significant systematic shift in the mag- 
nitudes was found, and the standard deviation of the retrieved 
magnitude distribution is adopted as the measured error at the 
respective output magnitude. 

On the basis of these simulations, we thus model the errors 
for the MSIS msoo? photometry, which increase approximately 
linearly towards fainter magnitudes, by 

e~0.25(m5oo7-28.5) [29.0,30.4]. (2) 

We then evaluate a completeness correction function, using the 
fraction of objects retrieved at each magnitude as these become 
fainter. This fraction is nearly 100% at 29.0 mag, the apparent 
magnitude of the PNLF bright cutoff at 5 1 Mpc, and decreases 
linearly down to 10-20% at 30.4 mag, the detection limit mag- 
nitude of our observations. We model this dependence by 

ri ifm5oo7 < 29.0, 

/, ~ <^ 0.64(-m5oo7 + 30.55) if 29.0 < m5oo7 < 30.4, (3) 
[0 ifm5oo7 > 30.4. 

The error distribution and the completeness function are used in 
Section Q below to perform simulations of the LOSVD for the 
PN sample. 
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Fig. 1. PNs in the Hydra I cluster core. Left panel: the PN velocity-magnitude distribution. The black crosses show the entire 
sample of 56 PN candidates. The blue and red lines are the measured transmission curves of the blue and the red filter, respectively, 
normalized so that the maximum transmission is near the theoretical bright cutoff of the PNLF at the distance of Hydra I. Central 
panel: the PN LOSVD (black histogram). The bins in velocity are 270 km s^^ wide. The blue and the red solid lines show again the 
suitably normalized transmission curves of the blue and red filters. The vertical magenta, green and gray lines in both panels mark the 
systemic velocity of Hydra I, NGC 3311 and NGC 3309, respectively. Right panel: Spatial distribution of the PNs (black diamonds) 
in the MSIS field. The field is centered on NGC 33 11 and has size ^ 100 x 100 kpc'^; north is up and east to the left. The two open 
triangles indicate the positions of NGC 331 1 (center) and NGC 3309 (upper right). The PN indicated by the gray symbol is the only 
object compatible with a PN bound to NGC 3309, based on its position on the sky and LOS velocity, v gj.^y pj^j = 4422kms^^. 



6. The PN sample in Hydra I 

Our PN catalog for the central (100 kpc)^ of the Hydra I 
cluster contains 56 candidates, for which we measure vloSi 
xpN, ypN and TO5007- The detected PN velocities cover a 
range from 970kms~^ to 6400 kms^^ with fluxes from 2.2 x 
10~^®ergcm~^s^^ to 7.6 x 10~^®ergcm^-^s^^. The detected 
sample of objects have a magnitude distribution compatible with 
the PNLF at the distance of Hydra I; see also Section lTTI 

The magnitude-velocity plane - The properties of the PN 
sample in the velocity-magnitude plane are shown in the left 
panel of Fig.[fl In this plot, the apparent magnitude of the PNLF 
bright cutoff at the distance of the Hydra I cluster corresponds to 
a horizontal line at 29.0 mag. The blue and red lines are the filter 
transmission curves, as measured from the spectra, normalized 
so that the maximum transmission occurs near the PNLF bright 
cutoff. The PNs are indeed all fainter than msooy = 29.0 and 
extend to the detection limit magnitude, m ^jj. This is slightly 
different for the two filters; the faintest PNs detected through the 
blue filter have m g ^jj = 30.45, and those detected with the red 
filter have m ^jj = 30.3. 

The PN LOSVD - The measuredLOSVD of the PN sample is 
shown by the black histogram in the central panel of Fig. [T] The 
velocity window covered by the two filters is also shown and the 
systemic velocities of Hydra I, NGC 3311 and NGC 3309 (see 
Section|2]l are indicated by the magenta, green and gray vertical 
lines, respectively. These velocities fall in the middle of the ve- 
locity window allowed by the filters, where both filters overlap. 
The mean velocity of all PN candidates is vpNs — 3840 km s^^ 
and the standard deviation is rmsp^rs = 1390 km s"^ The dis- 
tribution is highly non Gaussian and dominated by several in- 
dividual components. The main peak appears in the range of 
velocities from 2400 to 4400 km s^^ and its maximum is at 
~ 3100 kms^^, within Icrjjy of the systemic velocity of the 

^ Thi s plot is based on more accurate photometry than and updates 
Fig. 1 of iVentimiglia et alj d2008l) . 



Hydra I cluster In the blue filter velocity range there is a sec- 
ondary peak at ^ 1800 km s~^ that falls 2 — 3(Tjjy from the 
systemic velocity of Hydra I. This blue peak may contain a few 
background galaxy contaminants, as discussed in Section 15.11 
above. Finally a red peak at ~ 5000 km s^^ within ^ 2(Tj_[y of 
the cluster mean velocity is detected in the velocity interval from 
4600 to 5400 kms^^, and there are some PNs with even higher 
LOS velocities. 

The spatial distribution of the PNs - The locations of the 
detected PNs on the sky are shown in the right panel of Fig. [T] 
Their spatial distribution can be characterized as follows: 

- most PNs follow an elongated north-south distribution ap- 
proximately centered on NGC 3311; 

- there is no secondary high density concentration around 
NGC 3309. Only one PN, indicated by the gray symbol in 
the right plot of Fig. [1] has a combination of velocity and 
position that are compatible with a PN bound to the halo of 
NGC 3309; 

- a possibly separate concentration of PNs is present in the 
northeastern corner of the field. 

We summarize our main results so far: 

1 . The PN candidates detected in the MSIS field have luminosi- 
ties consistent with a population of PNs at the distance of the 
Hydra I cluster. 

2. The distribution of PNs in the MSIS field is centered on 
NGC 3311. Only one candidate is consistent with being 
bound to NGC 3309, even though NGC 3309 is of compa- 
rable luminosity to NGC 331 1 and, on account of the X-ray 
results (see Section |2]), is likely located in the inner parts of 
the cluster within the dense ICL, at similar distance from us 
as NGC 3311. 

3. There is no evidence of a single, well-mixed distribution of 
PNs in the central 100 kpc of the Hydra I cluster, contrary 
to what one would expect from the dynamically relaxed ap- 
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pearance of the X-ray emission. Instead, the observed PNs 
separate into three major velocity components. 

7. Kinematic substructures and a parameter for 
tlie observed PN sample in Hydra I: 
comparison with a simulated MSIS model 

At this point, we would like to reinforce the last point by com- 
paring the observed velocity distribution with a simple model. 
The model is obtained by assuming a phase-mixed PN popu- 
lation placed at the distance (51 Mpc) and mean recession ve- 
locity of NGC 3311, and simulating its line-of-sight velocity 
distribution by convolving with the MSIS instrumental set up. 
The velocity dispersion of the PN population is taken to be 
4 64kms~^, the highest val ue measured from the long-slit data 
in lVentimiglia et al.l (l2010bl) . In this way we can test more quan- 
titatively whether the observed multipeaked LOSVD for PNs in 
our field is biased by the MSIS observational set-up or whether it 
provides evidence of un-mixed components in the Hydra I clus- 
ter core. 



7. 1. Predicting tlie luminosity function and LOSVD witli 
MSIS for a model PN population 

We first characterize the model in terms of the intrinsic luminos- 
ity function and LOSVD of the PN population. Then we describe 
the steps required to predict the coiTesponding TO5007 magnitude 
vs. LOS velocity diagram and LOSVD that would be measured 
with the MSIS set up. In the next subsection we compare the 
results obtained with the observed Hydra I PN sample. 

Model for the intrinsic PN population - The intrinsic 
PNLF can be approxi mated by the analytical function given by 
ICiardulloetal]([l989h : 



N{m) = Ce 



0.307 m 



1 



p3(m*- 



(4) 



where m is the observed magnitude, m* — 29.0 is the apparent 
magnitude of the bright cutoff at the adopted distance of NGC 
3311, and C is a multiplicative factor. The integral of N{m) 
from m* to m* + 8 gives the total number of PN associated 
with the bolometric luminosity of the parent stellar population 
(NpN in Eq. [B, and the C par ameter can be relate d to the 
luminosity-specific PN number a dBuzzoni et al ] l2006h . For our 
model we distribute the magnitudes of a PN population accord- 
ing to a very similar form ula fitted by Mendez to the results of 
iMendez & Soffneii(ll997h . 

Next we assume that this PN population is dynamically 
phase-mixed and that its intrinsic LOSVD is given by a Gaussian 
centered on the systemic velocity of NGC 33 1 1 , v. 



G(v) 



1 



o" corev 27r 



exp 



(v-v) 



21 



2ai 



(5) 



where here we adopt v ~ 3830 km s ^ (IVentimiglia et al .1201 Obi 
corrected to the filter frame), and for the velocity dispersion we 
take cr core — 464kms~^, the highest value measured from 
the long-slit data in this paper This approximates the velocity 
dispersion for the intracluster component in the outer halo of 
NGC 3 3 1 1, at central distances --20-30 kpc ("Ventimiglia et al.' 

We will consider the magnitude-velocity diagram and 
the LOSVD as histograms in velocity; then in each velocity bin 
Av, , the number of PNs is 



where G(v) is normalized so that J^i Av^ — 1. 

Simulating the MSIS observations - The magnitude-velocity 
diagram for such a model population is modified by a number 
of effects in the MSIS observations, which we simulate as de- 
scribed below. The MSIS simulation procedure implements the 
following steps: 

- the through-slit convolution of the PNLF; 

- the convolution with the filter transmission; 

- the photometric error convolution; 

- the completeness correction; 

- the computation of the LOSVD. 

The through-slit PNLF - The MSIS technique is a blind sur- 
vey technique. Therefore the positions of the slits on the sky are 
not centered on the detected objects, and the further away an 
object is from the center of its slit, the fainter it becomes. This 
effect is a function of both seeing and slit width, and it modifies 
the functional form of the PNLF, which needs to be accounted 
for when using the LF from MSIS PN detected samples. 

In principle, some PNs may be detected in two adjacent slits 
of the mask, and this would need to be corrected for. However, 
at the depth of the present Hydra I survey this effect is not im- 
portant for the predicted PNLF, and indeed no such object has 
been found in the sample. 

Given a "true" PNLF LF(m), the ''through slit PNLF" 
sLF(m) can easily be computed, and depends on slit width and 
seeing; for further details see Gerhard et al. (2011, in prepa- 
ration). The effect of the through-slit correction is to shift the 
sLF(m) faintwards in the observable bright part, compared to 
the "true" PNLF. 

Convolution with filter transmission - When the filter trans- 
mission T(vi) is less then 1 (100%), it shifts the through-slit 
PNLF to fainter magnitudes. The Am depends on the value of 
the filter transmission curve at the wavelength A or equivalent 
binned velocity v^, and is equal to ATO(vi) = — 2.5 logT(vi). 
The resulting instrumental PNLF, the distribution of source 
magnitudes before detection, becomes velocity dependent, i.e., 
iLF(m,Vj). 

For the present MSIS Hydra I observations, the combined 
filter transmission curve from both filters is defined as 



T(v), = max[Tg(v,),TR(v,)] 



(7) 



LF{w,) ~ N{m) G(vO Av, 



(6) 



where B and R denote the blue and red filters. It is 1 where the 
transmission is 100%, approximately from ~ 1500 kms^^ to ^ 
3300 kms"^ and from - 4200 kms"^ to - 6300 kms"^ it is 
< 1 in the filter gap around ^ 3800 km s^^ and at the low and 
high velocity ends of the observed range. 

Photometric error convolution - Once the instrumental LF 
iLF(TO,Vi) is computed, it must be convolved with the photo- 
metric errors which, for the case of the Hydra I observations, 
are given by the linear function in Eq. |2] Because of the photo- 
metric eiTors, PNs that are intrinsically fainter than the detection 
limit (here mag ~ 30.4) may be detected if they happen to fall 
on a positive noise peak on the CCD image, and PNs that are 
intrinsically brighter than mag ^ 30.4 may be lost from the sam- 
ple. Generally, because the through-slit PNLF sLF(m) increases 
towards fainter magnitudes, the photometric errors scatter more 
faint objects to brighter magnitudes than vice-versa; so the effect 
of the convolution is to shift the PNLF to brighter magnitudes 
again. 

Completeness correction - The completeness correction at 
a given observed magnitude is a multiplicative function which 
accounts for the decreasing fraction of PNs at fainter magnitudes 
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detected against the noise on the MSIS image. For the case at 
hand it is given in Eq.[3] After the last two steps, we aiTive at the 
final "MSIS PNLF', MSLF(m) for short. 

Computation of the simulated LOSVD - For each velocity 
bin the MSLF(r7i, v^) is integrated between the apparent magni- 
tude of the PNLF bright cut off (m* = 29.0 for Hydra I) and 
the detection limit magnitude in the relevant filter, m ^ (see 
Section |6]l, to obtain the "observed" cumulative number of PNs 
in each velocity bin: 

/•nif.dl 

^MSIS(V«)=/ MSLF(m,v,)dm. (8) 

The most cumbersome step in this procedure is the cor- 
rection for the filter transmission, because it makes the final 
MSLF(TO,Vi) velocity-dependent. It must correctly be applied 
before the convolution with the photometric eiTors, because the 
latter depend on the flux measured at certain positions on the 
CCD. So the errors on the through-slit magnitudes depend on 
the filter transmission values of the PNs. 

However, we have found that the observed MSLF for the 
Hydra I PN sample, when obtained from wavelength regions 
where the filter transmission is ^ 100%, is very similar to the 
one obtained by summing over the entire filter bandpass. The ef- 
fect of the velocity dependence on the overall MSLF must there- 
fore be small, and for the comparison of simulated and measured 
LOSVDs below we have therefore applied the filter transmission 
only after the error convolution and completeness correction. 

Before we discuss the LOSVD obtained from the complete 
model, we show in Fig. |2]the predicted cumulative luminosity 
function resulting from error convolution, completeness correc- 
tion, and filter transmission correction of the through-slit lumi- 
nosity function, weighting by the number of observed PNs in 
each velocity bin. Also shown in Fig. |2]is the cumulative his- 
togram of the 7715007 magnitudes for the 56 observed PNs in the 
MSIS field. With a cutoff magnitude of 29.0 the model fits the 
observed histogram fairly well; however, this is not a formal best 
fit to the distance. The important point shown by Fig.|2]is that the 
observed MSIS luminosity function of the PN emission sources 
in the Hydra cluster core is evidently consistent with a popula- 
tion of PNs at ~ 50 Mpc distance. 




Fig. 2. Cumulative luminosity function predicted for the present 
MSIS observations and the nominal cutoff magnitude of the 
Hydra I cluster, 29.0 (full red line, see text), compared with the 
cumulative histogram of the observed msooy magnitudes. 



7.2. Reality of observed kinematic substructures 

The simulated MSIS LOSVD given by TVmsis (vi ) for the simple 
Gaussian velocity distribution model and luminosity function of 
Eq.|4]is shown as the green histogram in Fig. [3] with the observed 
PN LOSVD overplotted in black. The simulated MSIS LOSVD 
is scaled to approximately match the observed Hydra I sample 
in the central velocity bins. 




1.' ,1 ,|J U-i--l^ 

1000 2000 3000 4000 5000 6000 

Fig. 3. LOSVD for the Hydra I PN sample from Fig. □ (black 
histogram), compared with a simulated MSIS LOSVD (green 
histogram) for a Gaussian velocity distribution with (Tcoio ~ 
464 kms^^; see text for further details. The blue-red solid line 
shows the combined filter transmission curve as given in Eq. [T] 
The vertical magenta, green and gray lines mark the systemic 
velocity of Hydra I, NGC 33 1 1 and NGC 3309, respectively. 



The comparison between the simulated LOSVD and the 
Hydra I PN LOSVD in Fig.[3]identifies the central peak at about 
3100 km s"^ in the observed PN LOSVD with that of the PN 
population associated with the stellar halo around NGC 331 1 in 
the cluster core, with ctcoic ^ 500 kms^^. The mean Vcorc and 
Ccoro of this component are approximately consistent with those 
of the intracluster light hal o of NGC 3311 de r ived fro m the long- 
slit kinematic analysis in IVentimi^ha et alJ (l20I0bl) . However, 
the asymmetry and offset of the peak of the observed histogram 
(by several 100 km s^^) relative to the MSIS convolved model 
centered at the systemic velocity of NGC 3311 appear signifi- 
cant (CTcoro/ V -^coro — 100 km s~^), arguing for some real asym- 
metry of the central velocity component. We shall refer to the 
central peak in the Hydra I PN LOSVD in Fig. [3] as the central 
ICL component. 

Two additional velocity peaks are seen in the LOSVD in 
Fig-El one near 1800 kms~^ and one at ^ 5000 kms^^, which 
do not have any correspondence with the velocity distribution 
derived for the simulated MSIS model. These velocity compo- 
nents cannot be explained as artifacts of the MSIS set up, in 
particular, the filter gap in the Bh-R filter combination. We will 
refer to these two velocity components as secondary blue and 
red peaks, respectively. They reveal the possible presence of two 
kinematical substructures in the core of Abell 1060, whose ori- 
gins must be investigated further; see Section [8] 
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7.3. Low a-parameter in the core of Hydra I 

We now compare the number of observed PNs with the ex- 
pectations from the luminosity distribution and kinematics in 
and around NGC 3311. One issue is the absence of a clear 
subcomponent of PNs with velocity dispersion ^ 150 — 
250 kms^^, as would be expected from the central ^ 25" of 
NGC 331 1 (Ventimiglia et al. 2010b). It is known that FN sam- 
ples in elliptical galaxies are generally not complete in the cen- 
tral regions because of the increasing surface brightness pro- 
file; PNs are hard to detect against the image noise in the 
bright centers. E.g., in observations with the Planetary Nebula 
Spectrograph, the threshold surface brig htness is typically i n 
the range = 20 - 22 mag / arcsec^ (ICoccato et al .1120091) . 
In the current Hydra I data, the PN sample is severely incom- 
plete at fly = 21.0 mag/ arcsec^ (only two PNs are seen at 
/iy ^21.0 mag/ arcsec^, and si x at /xy > 2 1.5 mag / arcsec^). 
Referring to Fig. 13 of Mendez et al.l (1200 ih . we estimate that 
the current sample is not complete for /iy < 22.0 mag / arcsec^, 
which is reached at a distance of ~ 30" from the center of 
NGC 33 1 1 (Arnaboldi et al. 201 1, in preparation). At this radius, 
the projected velocity dispersion has risen to crA'33ii(30") ~ 
300 - 400 km s"^ (Ventimiglia etani20I0bh . Thus the PNs de- 
tected in this paper almost exclusively sample the hot (intraclus- 
ter) halo of NGC 3311. The cold inner galaxy component is not 
sampled. 

The second issue is the observed total number of PNs, given 
the detection limit, the instrumental set up and the light in NGC 
3311 and NGC 3309. Integrating the simulated MSIS luminos- 
ity function down to the detection limit of 30.4 mag, we obtain 
an effective a parameter for our observations of aMsis.Hy = 
0.012atot, where atot quantifies the total number of PNs 8 mag 
down the PNLlB This value is similar to q:o.5> the integrated 
value 0.5 mag down the PNLF. It is consistent with Fig.[T] even 
though in this figure PNs are seen up to 1.5 mag fainter than the 
nominal cutoff magnitude, because of (i) the shift towards fainter 
magnitudes due to the slit losses, and (ii) the completeness cor- 
rection (Eq.O. 

We can estimate the bolometric atot for NGC 3311 from 
its (FUV-V) color, th e relation between (F UV-V) and logai.o 
shown in Fig. 12 of ICoccato et al.l (20091) , and correcting to 
log atot by using Fig 8 of lBuzzoni etal.1 (l2006l) . The (FUV- 
V) color is determined from the Galex FUV magnitude and 
the V band magnitud e from RC3, b oth co rrected for extinc- 
tion, as described in ICoccato et alJ (l2009l Section 6.1). The 
resulting value, (FUV-V)=6.7, corresponds to logai.o = 1-1 
and log atot — —7.34. This is very similar to the value of 
log at. nt = —7.30 found for the Fornax cluster cD galaxy NGC 
1399 (iBuzzoni et al . 2006). Using the V band light profile of 
NGC 3311 measured in Arnaboldi et al. (2011, in preparation), 
and a bolometric correction of 0.85mag, we can then predict the 
expected cumulative number of PNs within radius R from the 
center of NGC 3311. This is shown as the red curve in Figure|4] 
after subtracting the luminosity within 20" which is not sampled 
by our MSIS observations. Also shown are the cumulative his- 
tograms of the observed number of PNs in the MSIS data, both 
for all PNs in the field, and for PNs with velocities in the central 
velocity component only. 

Fig.|4]shows that the total number of PNs detected in the field 
falls short of the number predicted from the luminosity profile by 
a factor ~ 4. Outside ^ 100" , the number of PNs with velocities 



consistent with the central ICL halo of NGC 33 1 1 is a factor ^ 2 
lower than the number of all PNs. Clearly therefore, some of the 
light at these radii is in a component different from the phase- 
mixed central ICL halo, but the amount is uncertain because we 
do not know whether the luminosity-specific a-parameter of this 
component is similarly low as for the NGC 331 1 ICL halo. For 
example, agreement between observed and predicted PN num- 
bers could be achieved by scaling only the NGC 3311 halo com- 
ponent by a factor ^ 6. On the other hand, scaling only an outer 
component will not work, because the discrepancy in Fig.|4]is 
already seen at small radii. Thus we can conclude that the a- 
parameter of the NGC 331 1 ICL halo is low by a factor 4—6. 

Such an anomalous specific PN number density requires an 
explanation. One possibility is that the stellar population in the 
halo of NGC 331 1 is unusually PN poor; this will need studying 
the stellar population in the galaxy outskirts. A second possibil- 
ity is that the ram pressure against the hot X-ray emitting gas 
in the halo of NGC 331 1 is high enough to severely shorten the 



m the halo of JNUC Hi I is high enough to severely shorten the 
lifetime of the PNs dPopita et all 120001; IVillaver & S tanghellinil 
'2005'). In their simulations. ^ Villaver & Stanghellinil (12005 ) con- 
sider a gaseous medium of density n = 10~* cm~^ and a rela- 
tive velocity of 1000 kms^^. They find that the inner PN shell 
is not significantly affected by the ram pressure stripping during 
the PN lifetime, and because the inner shell dominates the line 
emission in their model, the PN visibility lifetime is therefore not 
shortened relative to an undisturbed PN. However, with a density 
of the ICM inside 5' around NGC 3311 of - 6 x lO^^ cm-^, 
and a typical velocity of \/3 x 450 kms^^ ~ 800kms~^ the 
ram pressure on the NGC 3311 is ^40 times stronger than in 
their simulated case, so the ram pressure effects could be much 
stronger. Unfortunately, simulations of the evolution of PNs in 
such dense media are not yet available, to our knowledge. 
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This value includes the light between adjacent slits for the normal- 
ization. 
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Fig. 4. Observed and predicted cumulative PN numbers, as a 
function of radial distance from the center of NGC 3311. The 
green line shows the cumulative number of PNs associated with 
the central ICL halo of NGC 331 1, based on their velocities. The 
black line shows the cumulative number of all PNs, without ve- 
locity selection. The red curve shows the predicted cumulative 
number of PNs computed using the luminosity-specific parame- 
ter a estimated as explained in the text 

, the MSIS observational set-up, and the integrated bolometric 
luminosity in increasing circular apertures centered on NGC 331 1. 

If this explanation is correct, PNs should be most efficiently 
ram pressure stripped in the innermost, densest regions of the 
ICM. Hence in this case we would expect most of the observed 
PNs to be located in the outermost halo of NGC 3311, even 
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those projected onto the inner parts of our MSIS field. At these 
outer radii, dynamical time-scales are longer, and phase-mixing 
should be less complete. This would fit well with the unmixed 
kinematics and spatial distribution of the observed sample (see 
also next Section). 

The third issue is that we do not see a concentration of PNs 
around NGC 3309. As shown in Section |6] only one PN in the 
sample, shown by the gray symbol in the right panel of Fig. [T] 
has both position and LOS velocity compatible with being bound 
to NGC 3309. Whereas using the relative total luminosities of 
NGC 3309 and NGC 3311 to scale the number of PNs associ- 
ated with the main LOS velocity component for NGC 3311 in 
Fig.[3](i.e., 27 PNs), we would expect about 11 PNs associated 
with the light of NGC 3309 if both galaxies were at the same dis- 
tance. There are two possible explanations for this fact. One is 
that NGC 3309 is at significantly larger distance than NGC 3311, 
such that even PNs at the bright cutoff would be difficult to see. 
However, a simple calculation shows that then NGC 3309 would 
be put at ^ 70 Mpc well outside the cluster, at variance with X- 
ray observations finding that its gas atmosphere is confined by 
the ICM pressure (see Section|2]i. The second possibility is that, 
similarly as for NGC 3311, also the PNs in NGC 3309 may be 
severely ram pressure stripped by the galaxy's motion through 
the dense ICM in the cluster core. This would require that NGC 
3309 moves rapidly through the cluster core, and is physically 
rather close to NGC 3311. Again simulations would be needed 
to check this quantitatively. 

8. The substructures in the Hydra I cluster core 

We now turn to a more general discussion of the spatial distri- 
bution and kinematics of PNs and galaxies in the central region 
of the cluster ICL is believed to originate from galaxies, so it 
is interesting to ask whether the phase-space substructures seen 
in the distribution of the PNs that trace the ICL has some cor- 
respondence to similar structures in the distribution of cluster 
galaxies. Thus we want to investigate the spatial distributions of 
the PNs associated with the velocity subcomponents in the PN 
LOSVD discussed earlier, and compare them with the spatial 
distribution of Hydra I galaxies in similar velocity intervals. In 
this way, we may obtain a better understanding of the dynamical 
evolution of the galaxies in the cluster core, and of the relevance 
of cluster substructures for the origin of the diffuse cluster light. 

8. 1. Spatial distributions of tlie PN velocity components 

We first consider the spatial distribution of the PNs associated 
with the different velocity components in the PN LOSVD. This 
is shown in the three panels of Fig. |5] divided according to the 
classification in Sect. 17.21 Each panel covers a region of 6.8 x 
6.8arcmin^ ~ 100 x 100 kpc^ centered on NGC 3311. 

PNs associated with the central ICL component (middle 
panel of Fig.|5]l can be divided into two spatial structures. There 
is a prominent PN group concentrated, as expected, around 
NGC 331 1, and an elongated east-west distribution in the north- 
ern part of the FoV. By contrast, we see a low PN density region 
in the southern part of the MSIS field. 

Such a north/south asymmetry is seen also in the spatial dis- 
tribution of the galaxies. Fig. |6] displays a larger area, 20 x 
20 arcmin^, which includes the MSIS field studied in this work, 
as indicated by the orange square. We can see from the two pan- 
els (photo, and schematic) that NGC 331 1 and NGC 3309 dom- 
inate the center of the MSIS field, that there is a high density of 



bright galaxies in the northern part of the field, but a deficit of 
galaxies to the south of NGC 3311. 

The spatial distribution of the PNs associated with the sec- 
ondary red peak in the PN LOSVD is shown in the right panel of 
Fig.|5] It has a north/south elongation, apparently extending fur- 
ther towards the south of NGC 331 1 than the central ICL com- 
ponent, with a high density region north/east of NGC 331 1. 

Finally, the spatial distribution of the PNs associated with the 
secondary blue component at 1800 kms^^ (left panel ofFig.|5]) 
also appears elongated along the north/south direction, but the 
smaller number of objects in this subsample makes inferring 
their spatial structure more difficult. 

In summary, there is little evidence of a spherically symmet- 
ric well-mixed distribution of PNs in the outer halo of NGC 3311 
in the cluster core. Several velocity components are seen, and 
even the central ICL component centered on NGC 3311 shows 
signs of spatial substructures. 

8.2. Spatial and velocity distribution of Hydra I galaxies: 
comparison witfi the PNs sample 

The spatial distri b ution of the g a laxies from 
IChristlein & Zabludof j (l2003h : iMiseeld et all (l2008h in the 
central 20 x 20arcmin^ centered on NGC 3311 is shown in 
Fig. |6] We would like to analyze their phase-space distribution 
by dividing into the same velocity components as identified in 
the PN LOSVD. Therefore, in the image on the left the bright 
galaxies are encircled with the colors of the PN components in 
Fig. |5] and in the right panel all galaxies in the field are shown 
schematically as squares and crosses with the same color code 
for these velocity bins. NGC 3311 and NGC 3309 are marked 
in the center of the MSIS field (orange square). 

In Fig. |2l the left panel shows the velocity distribution of 
all the galaxies in the 20 x 20 arcmin^ region centered on 
NGC 3311. In the right panel, the velocity histograms for the 
bright galaxies (m^ < 15.37, violet color) and dwarf galaxies 
(m/f > 15.37, green color) are shown separately. 

The LOSVD for the Hydra I galaxies covers the same ve- 
locities as for the PN sample. If we select only galaxies in 
the range of velocities of the PNs in the central ICL compo- 
nent, from 2800 kms"^ to 4450kms"\ their LOSVD is con- 
sistent with a Gaussian distribution centered at a velocity of 
3723 ± 100 kms"^ with a dispersion of 542 ± 80kms~\ This 
confirms re sults from long-slit kine matics in the outer halo of 
NGC 331 1 (IVentimiglia et al.ll2010bh . where the velocity disper- 
sion was found to increase to ~ 465 km s~^ at ^ 70" radius, 64% 
of the velocity dispersion of all cluster galaxies. 

This subsample of galaxies also has an interesting spa- 
tial distribution: the central 6.8 x 6.8 arcmin^ region of the 
cluster (the MSIS field), while dominated by NGC 3311 and 
NGC 3309, contains no other Hydra I galaxies with these ve- 
locities. Whereas outside this region, they appear uniformly dis- 
tributed over the field (see the green squares and crosses in the 
right panel of Fig.|6ll. NGC 331 1 is at the center of the distribu- 
tion of these galaxies both in space and in velocity. The distribu- 
tion of these galaxies, as well as the similarity of their velocity 
dispersion with that measured in the ha l o of NG C 3311, supports 
the interpretation of IVentimiglia et al] (I2010bh that the halo of 
NGC 331 1 is dominated by intracluster stars that have been torn 
from galaxies disrupted in the cluster core: galaxies that passed 
through the central 100 kpc of the cluster core at modest veloci- 
ties have all been disrupted. 
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Fig. 5. Left panel: Spatial distribution of the PNs associated with the blue secondary peak in the PN LOSVD (< 2800 km s~^). 
Central panel: Spatial distribution of the PNs associated with the central ICL component ( 2800 kms^^ to 4450 kms^^). Right 
panel: Spatial distribution of the PNs associated with the secondary red peak at > 4450 km s^^ in the PN LOSVD. The black 
triangles indicate NGC 33 11 (center) and NGC 3309 (north-west of center), respectively. North is up and east is to the left. 
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Fig. 6. Left panel: 20 x 20arcmin DSS image of the Hydra 1 cluster. The two bright galaxies at the field center are NGC 3311 
(center) and NGC 3309 (north-west of center). The blue circles indicate galaxies with Vgys < 2800 kms^^, the green circles 



galaxies with 2800 km s 



< Vs 



< 4450 km s ^ (only those within 10 arcmin around NGC 331 1 and with m/? > 15.37), and the 



red circles galaxies with Vgys > 4450 km s~ . Right panel: Spatial di stribution of Hy dra 1 galaxies in the same area of 20 arcmin 
cente red on NGC 3311. Sq uares indicate galaxies from the catalog of IChristlein & Z abludoff (2003) and crosses indicate galaxies 
from lMisgeld et all (l2008l) . The color of the symbols refers to the velocity components in the PN LOSVD as described in Fig.|5] 
The two diamonds locate NGC 331 1 and NGC 3309. The orange square shows the FoV used in the FORS2 MSIS observations. 



By contrast, the galaxies with LOS velocities > 4450 km s~ 
as in the secondary red peak of the PN LOSVD are mostly lo- 
cated within the central 100 x 100 kpc^ region of the cluster 
(red squares and crosses in the right panel of Fig.|6]l. In this sub- 
sample, there are 14 galaxies in total, 5 are classified as bright 
galaxies and 9 are dwarfs, and 3 bright galaxies and 6 dwarfs fall 
within the MSIS FORS2 field. These 6 dwarfs are concentrated 
in the northeastern part of the halo of NGC 3311, in the same 



region occupied by many PNs associated with the secondary red 
peak. 

Finally, in this region there are only a few galaxies with a 
LOS velocity lower than 2800 km s^^, compatible with the sec- 
ondary blue peak in the PNs. They are 8 in total (blue squares 
and crosses in the right panel of Fig.|6]l. Only one of these falls 
on the boundary of the central 100 x 100 kpc^ region around 
NGC 331 1. One of these galaxies is the giant spiral NGC 3312, 
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Fig. 7. Left panel: Histogram showing the LOS VD of all galaxies from the catalog of ' Christlein & Zabludoffl (12003 *) within an area 
of 20 arcmin in size centered around NGC 33 11. Right panel: The purple histogram indicates t he LOSVD for the br ight galaxies 
in this field, and the light green histogram the LOSVD for all dwarf galaxies from the catalog of iMisgeld et alJ (l2008h in the same 
area. 



south-east of NGC 3311. The others, including the spiral galaxy 
NGC 3314, are located at larger distances from NGC 331 1. 

8.3. Galaxy evolution and presence of substructures in 
the core of the Hydra I cluster 

The distribution of galaxy properties in clusters holds important 
information on galaxy evolution and the growth of galaxy clus- 
ters. In Sect. 18.21 we have discovered an apparent lack of gal- 
axies in the central 100 x 100 kpc^ region of the cluster core with 
velocities in the same range as covered by the cD h alo. A simi- 
lar re sult has been found in the NGC 5044 group ( Mendel et al.l 
|2009|) . A possible explanation is the tidal disruption of galaxies 
at small cluster-centric radii. Galaxies with LOS velocities in the 
range of the central ICL component of the Hydra I cluster are 
no longer seen in the central region of the cluster, because they 
were all disrupted in the past during close encounters with the 
lumino us galaxy and the dark matter distribution at the cluster 
center (' Faltenbacher & Mathewsll2005l) . Their former stars now 
contributes to the diffuse stellar component in the Hydra 1 core. 

Differently from the NGC 5044 group, however, we have 
found a number of dwarf galaxies with high velocities in the 
Hydra 1 core, with small cluster-centric radii (< 100 kpc). These 
dwarfs have LOS velocities higher than 4400 km s^^ and seem 
to form a well defined substructure both in velocity and spatial 
distribution. We speculate that these galaxies are falling through 
the cluster core and are not yet disrupted by the tidal interaction 
with NGC 33 1 1 . The PNLF of the subsample of PNs associated 
spatially and in velocity with this substructure places it almost 
at the same distance as the central ICL component of the cluster: 
this group of galaxies may indeed now be on the point of close 
encounter with NGC 331 1 in the cluster center. 

Finally, we have found a correlation between the PNs con- 
tributing to the secondary blue peak of PN LOSVD, and 8 
galaxies with a LOS velocity lower than 2800 km s^^. Among 
these galaxies is the large spiral galaxy NGC 3312 (Vs^^ = 
2761 km s~^) as well as NGC 3314 (v^„, = 2 795km s~^). 
iFitchett & Merritj (Il988l) and 'McMaho n et"an (1 19921) have 
claimed the presence of a foreground group of galaxies asso- 
ciated with these spirals. Unfortunately due to the small area 
covered by the current MSIS survey, it is difficult to determine 



unambiguously whether the low velocity PNs (which from their 
PNLF are at the distance of the cluster) are associated with these 
galaxies. A PN survey covering the region between NGC 331 1 
and NGC 3312 may provide a definite answer to this question. 

9. Summary and Conclusions 

Using multi-slit imaging spectroscopy with FORS2 on VLT- 
UTl, we have studied a sample of 56 planetary nebula (PN) 
candidates in the Hydra I cluster at 50Mpc distance, target- 
ing a region 100 x 100 kpc^ centered on the cluster cD galaxy, 
NGC 331 1. The MSIS technique allows us detect these emission 
sources and measure their velocities, positions and magnitudes 
with a single observation. 

PN candidates are defined as unresolved emission sources 
without measurable continuum. Emission sources that are either 
resolved spatially or in wavelength or have a detected contin- 
uum are classified as background galaxies; see Ventimigl ia et al.l 
(1201 Oal) . We show that the luminosity function of the PN can- 
didates is as expected for a population of PNs at the distance 
of Hydra I. Moreover, almost all the detected background gal- 
axies occur in the velocity range between 1000 kms^^ and 
2800kms"\ blue-shifted by > 900 kms"^ with respect to the 
mean recession velocity of the Hydra I cluster. From these facts 
we conclude that any residual contamination of the PN sample 
by background galaxies with undetectable continuum must be 
small and restricted to the velocity range given. 

The luminosity-specific number density a inferred from the 
PN sample and the luminosity of diffuse light around NGC 3311 
is a factor ^ 4 lower than expected, even if we compare with 
the low a value determined from the (FUV-V) color which is 
one of the lowest for elliptical galaxies. A possible interpretation 
is that ram pressure stripping by the dense, hot X-ray emitting 
intracluster medium in the center of the cluster core around NGC 
3311 dramatically shortens the lifetime of the PN phase. This 
also seems the most likely explanation for the observed lack of 
PNs bound to NGC 3309, the other giant elliptical galaxy in the 
Hydra I core. 

The line-of-sight velocity distribution (LOSVD) of the ob- 
served PNs shows at least three separate peaks, and their phase- 
space distribution is inconsistent with a single well-mixed intra- 
cluster distribution. One peak, which we term the central intra- 
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cluster component, is broadly consistent with the outward con- 
tinuation of the intracluster halo of NGC 3311, which was earlier 
shown to have a velocity disper sion of ^ 470 km s^^ at radii of 
> 50" (IVentimiglia et al.l2010bl) . Simulating MSIS observations 
for a Gaussian intrinsic LOSVD with ^470 km s^^ centered on 
the systemic velocity of NGC 3311 has additionally shown sig- 
nificant residual asymmetries, suggesting that also this central 
component is not completely phase-mixed in the central cluster 
potential. 

Many cluster galaxies are found in the LOS velocity range 
associated with this central intracluster component (2800 km s^ ^ 
to 4450 km s~^), but none in the central 100 x 100 kpc^ around 
NGC 331 1. We suggest that the missing galaxies have been dis- 
rupted by the gravitational field of NGC 331 1 and the surround- 
ing cluster dark matter, and that their light has been added to 
the diffuse intracluster halo of NGC 3311 which is traced by the 
PNs. 

The second main peak in the PN LOSVD is centered at 
5000 km s^^, some 1200 km s^^ to the red of the main compo- 
nent. In the same velocity range, a number of dwarf galaxies 
are seen, which are projected onto the central 100 x 100 kpc^ 
around NGC 3311 where also the PNs are located. We suggest 
that the PNs and the galaxies in this red peak of the LOSVD are 
linked, i.e., on similar orbits through the cluster core, indicat- 
ing that the galaxies have been partially disrupted and the tidal 
debris is traced by the PNs. This will be the subject of a fur- 
ther study based on deep photometry (Arnaboldi et al. 2011, in 
preparation). 

Finally, a third, blue peak in the PN LOSVD is seen at 
~ 1800 kms~^. The spatial distribution of these PNs is elon- 
gated in the same sense as for the other two components in the 
cluster core, but the number of sources with these velocities is 
smaller and a few of them might be unresolved background gal- 
axies. This makes it difficult to establish a robust association be- 
tween these PNs and cluster galaxies, such as the group related 
to the spiral NGC 3312. A larger survey area would be needed 
to establish such a link. 

In summary, from this study of the kinematics of diffuse light 
in the Hydra 1 cluster core with PNs, and the comparison with 
the projected phase-space distribution of galaxies, we infer that: 
(1) The intracluster stellar population in the Hydra 1 cluster is 
not well-mixed, even though this cluster is believed to be the 
prototype of an evolved and dynamically relaxed cluster based 
on X-ray indicators. (2) The build-up of diffuse intracluster light 
and of the cD halo of NGC 3311 are ongoing, through the ac- 
cretion of material from galaxies falling into the cluster core and 
tidally interacting with its potential well. 

Acknowledgements. The authors thank the ESO VLT staff for their support dur- 
ing the MSIS observations. They also thank L. Coccato, K.C. Freeman and 
E. lodice for useful discussions. This research has made use of the Gemini data 
archive and the NASA/IPAC Extragalactic Database (NED) operated by the Jet 
Propulsion Laboratory, California Institute of Technology. 



References 

Adami, C, Slezak, E., Dun-et, F., et al. 2005, A&A, 429, 39 
Aguerri, J. A. L., Gerhard, O. E., Arnaboldi, M., et al. 2005, AJ, 129, 2585 
Ai-naboldi, M., Doherty, M., Gerhard, O., et al. 2008, ApJ, 674, L17 
Arnaboldi, M. & Gerhard, O. 2010, Ai'Xiv e-prints, 1001.3523 
Arnaboldi, M., Gerhard, O., Aguerri, J. A. L., et al. 2004, ApJ, 614, L33 
Arnaboldi, M., Gerhard, O., Okamura, S., et al. 2007, PASJ, 59, 419 
Buzzoni, A., Arnaboldi, M., & Corradi, R. L. M. 2006, MNRAS, 368, 877 
Castro-Rodriguez, N., Arnaboldi, M., Aguerri, J. A. L., et al. 2009, A&A, 507, 
621 

Christlein, D. & Zabludoff, A. 1. 2003, ApJ, 591, 764 



Ciardullo, R., Jacoby, G. H., Feldmeier, J. J., & Bardett, R. E. 1998, ApJ, 492, 
62 

Ciardullo, R., Jacoby, G. H., Ford, H. C, & Neill, J. D. 1989, ApJ, 339, 53 
Ciardullo, R., Sigurdsson, S., Feldmeier, J. J., & Jacoby, G. H. 2005, ApJ, 629, 
499 

Coccato, L., Arnaboldi, M., Gerhard, O., et al. 2010, A&A, 519, A95+ 

Coccato, L., Gerhard, O., Arnaboldi, M., et al. 2009, MNRAS, 394, 1249 

Doherty, M., Arnaboldi, M., Das, P, et al. 2009, A&A, 502, 771 

Dolag, K., Murante, G., & Borgani, S. 2010, MNRAS, 405, 1544 

Dopita, M. A., Massaglia, S., Bodo, G., Arnaboldi, M., & Merluzzi, P. 

2000, in Astronomical Society of the Pacific Conference Series, Vol. 199, 

Asymmetiical Planetary Nebulae 11: From Origins to Microstructures, ed. 

J. H. Kastner, N. Soker, & S. Rappaport, 423-+ 
Durrell, P. R., Ciardullo, R., Feldmeier, J. J., Jacoby, G. H., & Sigurdsson, S. 

2002, ApJ, 570, 119 
Faltenbacher, A. & Mathews, W. G. 2005, MNRAS, 362, 498 
Feldmeier, J. J., Ciardullo, R., Jacoby, G. H., & Durrell, P R. 2004, ApJ, 615, 

196 

Fitchett, M. & Merritt, D. 1988, ApJ, 335, 18 

Furusho, T., Yamasaki, N. Y, Ohashi, T. S. R., et al. 2001, PASJ, 53, 421 
Gal- Yam, A., Maoz, D., Guhathakurta, P, & Filippenko, A. V. 2003, AJ, 125, 
1087 

Gerhard, O., Arnaboldi, M., Freeman, K. C, et al. 2005, ApJ, 621, L93 
Gerhard, O., Arnaboldi, M., Freeman, K. C, et al. 2007, A&A, 468, 815 
Gnedin, O. Y. 2003, ApJ, 589, 752 

Gonzalez, A. H., Zabludoff, A. I., & Zaritsky, D. 2005, ApJ, 618, 195 
Hausman, M. A. & Ostriker, J. P 1978, ApJ, 224, 320 

Hayakawa, A., Furusho, T, Yamasaki, N. Y., Ishida, M., & Ohashi, T. 2004, 
PASJ, 56, 743 

Hayakawa, A., Hoshino, A., Ishida, M., et al. 2006, PASJ, 58, 695 
Henize, K. G. & Westerlund, B. E. 1963, ApJ, 137, 747 

Herrmann, K. A., Ciardullo, R., Feldmeier, J. J., & Vinciguerra, M. 2008, ApJ, 
683, 630 

Hilker, M. 2002, in lAU Symposium, Vol. 207, Extragalactic Star Clusters, ed. 

D. P Geisler, E. K. Grebel, & D. Minniti, 281-+ 
Hilker, M. 2003, in Extragalactic Globular Cluster Systems, ed. M. Kissler-Patig, 

173-+ 

Jacoby, G. H. 1980, ApJS, 42, 1 
Jacoby, G. H. 1989, ApJ, 339, 39 

Kapferer, W., Schindler, S., Knollmann, S. R., & van Kampen, E. 2010, A&A, 
516, A41+ 

Kelson, D. D., Zabludoff, A. I., WilHams, K. A., et al. 2002, ApJ, 576, 720 
Krick, J. E. & Bernstein, R. A. 2007, AJ, 134, 466 
Lee, M. G., Park, H. S., & Hwang, H. S. 2010, Science, 328, 334 
McMahon, P M., van Gorkom, J. H., Richter, O., & Ferguson, H. C. 1992, AJ, 
103, 399 

McNeil, E. K., Arnaboldi, M., Freeman, K. C, et al. 2010, A&A, 518, A44+ 
Mendel, J. T, Proctor, R. N., Rasmussen, J., Brough, S., & Forbes, D. A. 2009, 

MNRAS, 396, 2103 
Mendez, R. H., Riffeser, A., Kudritzki, R.-P, et al. 2001, ApJ, 563, 135 
Mendez, R. H. & Soffner, T. 1997, A&A, 321, 898 
Mieske, S., Hilker, M., & Infante, L. 2005, A&A, 438, 103 
Mihos, J. C, Harding, P, Feldmeier, J., & Morrison, H. 2005, ApJ, 631, L41 
Misgeld, I., Mieske, S., & Hilker, M. 2008, A&A, 486, 697 
Moore, B., Katz, N., Lake, G., Dressier, A., & Oemler, A. 1996, Nature, 379, 

613 

Murante, G., Giovalli, M., Gerhard, O., et al. 2007, MNRAS, 377, 2 

Neill, J. D., Shara, M. M., & Oegerle, W. R. 2005, ApJ, 618, 692 

Pierini, D., Zibetti, S., Braglia, F, et al. 2008, A&A, 483, 727 

Puchwein, E., Springel, V, Sijacki, D., & Dolag, K. 2010, MNRAS, 406, 936 

Renzini, A. & Buzzoni, A. 1 986, in Astrophysics and Space Science Library, Vol. 

122, Spectral Evolution of Galaxies, ed. C. Chiosi & A. Renzini, 195-231 
Richstone, D. O. 1976, ApJ, 204, 642 

Rudick, C. S., Mihos, J. C, Frey, L. H., & McBride, C. K. 2009, ApJ, 699, 1518 
Rudick, C. S., Mihos, J. C, & McBride, C. 2006, ApJ, 648, 936 
Sato, K., Yamasaki, N. Y, Ishida, M., et al. 2007, PASJ, 59, 299 
Tamura, T, Makishima, K., Fukazawa, Y, Ikebe, Y, & Xu, H. 2000, ApJ, 535, 
602 

Vasterberg, A. R., Lindblad, P O., & Jorsater, S. 1991, A&A, 247, 335 
Ventimiglia, G., Arnaboldi, M., & Gerhard, O. 2008, Astronomische 

Nachrichten, 329, 1057 
Ventimiglia, G., Arnaboldi, M., & Gerhard, O. 2010a, La Rivista del Nuovo 

Cimento, in press 

Ventimiglia, G., Gerhard, O., Arnaboldi, M., & Coccato, L. 2010b, A&A, 520, 
L9+ 

Villaver, E. & Stanghellini, L. 2005, ApJ, 632, 854 

Williams, B. F, Ciardullo, R., Durrell, P R., et al. 2007, ApJ, 656, 756 

Willman, B., Governato, R, Wadsley, J., & Quinn, T. 2004, MNRAS, 355, 159 



12 



Ventimiglia et al.: Kinematics and origins of diffuse light in Hydra I cluster core 

Yamasaki, N. Y, Ohashi, T., & Furusho, T. 2002, ApJ, 578, 833 
Zibetti, S., White, S. D. M., Schneider, D. P., & Brinkmann, J. 2005, MNRAS, 
358, 949 



13 



